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We consider the running coupling from the four-gluon vertex in Landau gauge, SU(A'^c) Yang- 
Mills theory as given by a combination of dressing functions of the vertex and the gluon propagator. 
We determine these functions numerically from a coupled set of Dyson- Schwinger equations. We 
reproduce asymptotic freedom in the ultraviolet momentum region and find a coupling of order one 
at mid-momenta. In the infrared we find a nontrivial (i.e. nonzero) fixed point which is three orders 
of magnitude smaller than the corresponding fixed point in the coupling of the ghost-gluon vertex. 
This result explains why the Dyson- Schwinger and the functional renormalization group equations 
for the two point functions can agree in the infrared, although their structure is quite different. Our 
findings also support Zwanziger's notion of an infrared effective theory driven by the Faddeev-Popov 
determinant. 
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I. INTRODUCTION 



In recent years the running couphng of Yang-Mills theory has been investigated in a number of approaches; for 
a review see [T^. These include lattice QCD % i, i, i, S 0, 

analytic perturbation theory [1, the functional 
renormalization group [Tol. [Til [l^. Dyson-Schwinger equations [Tsl. [l^ fig and phenomcnological extractions from 
experiment [l^ • The goal of these investigations is an extension of our knowledge of the coupling from the large 
momentum region towards small momenta of the order of Aqcd and smaller. Perturbation theory alone, plagued by 
the problem of the Landau pole, is clearly insufhcient for this task. In this respect it seems remarkable that the mere 
improvement of the perturbation series by analyticity constraints leads to a well defined running coupling that freezes 
out in the infrared; see @ for a review of analytic perturbation theory. 

Infrared fixed points of the couplings of Yang-Mills theory have also been found in two functional approaches to 
QCD, the functional (or 'exact') renormalization group (FRG) and the framework of Dyson-Schwinger equations 
(DSEs); see [il[il,[2^ for reviews. In these approaches nonperturbative running couplings can be defined in terms 
of (gauge dependent) dressing functions of propagators and dressing functions of the primitively divergent vertices of 
the theory. The resulting expressions are renormalization group invariants but may be scheme dependent. In Landau 
gauge, the couplings from the ghost-gluon vertex, a^^"^'- , the three-gluon vertex, a^^ , and the four-gluon vertex, a^^, 
are given by [l5[ : 

„2 

gh-gli 



a^'^-y\p^) = 3^G'{p^)Z{p'), (1) 
47r 

'ip') - ^[r'^{p')rz'{p'), (2) 

'ip') = ^[r''ip')]z\p'). (3) 

Here g^/47r is the coupling at the renormalization point /i^, whereas Z{p^) denotes the dressing function of the gluon 
propagator and G{p^) the dressing of the ghost propagator D'^, i.e. 

= o^Ap'}={>^..-'-f)^. (4) 

The functions F^^ and F^f describe the nonperturbative dressing of the tree-level tensor structures of the three- and 
four-gluon vertices. The multiplicities of the various dressing functions in (IT])-([3]) are related to the number of legs 
of the corresponding vertex^. The three definitions of the coupling given in eqs. ([TJS]) correspond to three different 



Note that the ghost-gluon vertex is finite in Landau-gauge, which explains the absence of a corresponding dressing function in eq.JTJ. 
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renormalisation schemes. The resultin g co uphngs are related to each other by scale transformations and Slavnonv- 
Taylor identities as detailed e.g. in ref. [2l|. In this work we focus on a calculation of a^3{p^) and compare the result 
with the previously determined coupling a^^~f'(p^) [3, H^- 

One of the basic ingredients to the running coupling a'^^{p^) is the dressing function F^^ of the nonpcrturbative 
four-gluon vertex. An evaluation of this dressing together with a corresponding evaluation of the giuon propagator 
therefore allows to study the running of the coupling with momentum. However, there are also other reasons why the 
nonpcrturbative four-gluon vertex is an interesting object. First of all, this vertex is the only primitively divergent one 
that allows for the formation of bound state (glueball-) poles, a phenomenon usually restricted to higher, superficially 
convergent vertices. Second, the vertex describes quantum corrections to elementary gluon-gluon scattering, which 
might be important e.g. for the description of gluon-gluon interactions in the high temperature quark-gluon plasma 
phase of QCD. Third, a number of studies indicate [ij, [3, [H, US [13] that the infrared structure of the correlation 
functions of Yang-Mills theory is connected to the confining properties of the theory via the so called Gribov-Zwanziger 
scenario. Here, long ranged correlations are induced in the gauge fixed theory by effects from the first Gribov horizon 
in gauge field configuration space [131 ■ As we will see in the course of this work, the infrared behavior of the four-gluon 
vertex and the related running coupling provide additional support of this picture. 

This four-gluon vertex is a highly complex object due to its rich tensor structure generated by the four Lorentz 
and four color indices. As a consequence, this correlation function is very poorly understood so far. Lattice cal- 
culations of many-gluon Green's functions suffer from problems with statistics and consequently no definite results 
have been obtained so far. Within the functional continuum approach to Yang-Mills theory, early investigations of 
the vertex concentrated on the structure of its Dyson-Schwinger equations (see e.g. [2^), without aiming at actual 
solutions. Results on the one-loop level have been given e.g. in [2l|, [2^, [23|. An attempt to solve the vertex-DSE 
non-perturbatively has been made in [2^ [2^, [sO] within a selfconsistent expansion scheme in terms of couplings and 
power laws of momenta. 

In this work we are going beyond these results by a combination of analytical and numerical methods that allow to 
extract the dressing functions of the vertex without any prejudice to their functional form. In section |TT] we construct 
an approximation to the full Dyson-Schwinger equation of the vertex which reproduces the correct asymptotic behavior 
of the vertex as known from perturbation theory and infrared power counting methods . In section IIIII We give 
analytical expressions for the vertex in these two limits and discuss numerical results for all momenta in section IIVI 
For the running coupling Q!4g(p^) we find an infrared fixed point, which we discuss in section |Vl We explain why the 
smallness of this fixed point matches with results from the functional renormalization group and the notion of ghost 
dominance in the infrared. A summary and outlook concludes the paper. 



II. THE FOUR-GLUON VERTEX AND ITS DYSON-SCHWINGER EQUATION 
A. Nonperturbative structure of the four-gluon vertex 

As already mentioned above, the four-gluon vertex is a highly complicated object with four Lorentz- and four color 
indices. This complexity forces a two step procedure: one first works with a restricted subset of possible combinations 
of Lorentz- and color tensors. This reduced complexity allows for a first study of the most important properties of 
the vertex and its Dyson-Schwinger equation. On the basis of these results one can then attack the full problem in a 
second step. While we report on the first part of this program in this work, the second part is left for future studies. 
Of course, the success of such a procedure greatly depends on the choice of the restricted subset. A suitable selection 
has been suggested in [30i] and shall also be used here. 

The building blocks of the reduced tensor-structure are three Lorentz- and five color-tensors: 



^ibcd — ^abScd, G^^bcd — ^ac^bd, Cj^^bcd — ^ad^bc, 

^ abed ~ fabnfcdn^ ^ abed — facn fbdn- (6) 

This is the minimal subset of all possible tensor-structures, which has the following properties [soj : 



Furthermore, the bare four-gluon vertex is proportional to instead of g which leads to factors of F^^Z^ instead of the naive expectation 
Jr4gj2^4 fi-Qm ^^jjg number of legs in the vertex. 
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• It is dynamically closed under DSE and Bethe-Salpeter iterations, provided the only color-structure appearing 
in the three gluon vertex is fabc- 

• It closes under crossing operations. 

• It contains the structure of the bare four-gluon vertex. 

The last property of this subset allows for the representation of the high momentum limit of the vertex in this basis 
and also allows for the calculation of the relevant dressing function for the running coupling, eq. ([3]). 

From these tensors a basis of the linear space of Lorentz/color-tensors is constructed as a direct product 

^ {i,j)-abcd ~ ^abcd^{]) ' \' I 

where we abbreviate the various combinations as follows: 

Bi3 = L(3)C(3), Bi4 = i(3)C(4), Bi5 = L(3)C(5). 

Here the Lorentz- and color- indices are left implicit. The four-gluon vertex is then represented by: 

5 3 

''TafccT iPl,P2,P3) ^^Yl (Pl'P2,P3) T^!^jy,abcd^ (9) 

where the T'^^^'^ , , are elements of an orthonormal basis constructed from the elements i?i i r such that the tree- level 

{i,j):abcd 

vertex is included. The tensors of this basis can be found in appendix |Al The algebraic manipulations involved in 
the construction of this basis and also the one below have been performed with the use of FORM [3lj . 

Of course, the dressing functions Fy {pi,P2,P3) of the basis © are not completely independent. Bose symmetry of 
the four external vertex legs dictates interrelations between combinations of the F^ ipi,P2,P3)- This symmetry is of 
course reproduced by the exact vertex-DSE, although it is far from trivial how this works in detail, since one external 
leg is always connected with a bare internal vertex while the others are connected with dressed Green's functions. 
Thus any approximation to the full system is endangered to generate unsymmetric terms. These can (partly) be 
projected out by contraction with a reduced basis of tensor structures, which only include Bose symmetric objects. 
The construction of this reduced basis is described in appendix[Bl Here we only give the result in terms of the building 
blocks, eqs. ®,® ,®: 

= 108N^{N^-1) ( - ^4 + 2i?5 + 2B, Bio - B,, B,,) (10a) 

1/24 4 

Vo — i -Bi H Ba — -\- Bj — Bot 

_ 1207V2 + 72 V iNc 3iVc iNc 

2 2 2 \ 

H B10 + B13-I Bi4H Bi5 (10b) 

37Vc 3A^c 3iVc / 

1 



216 {N^ " + 7V2 + 4) 



iV2 -I- 6 2(7V2 + 1) 4(7V2 - 1) + 6 



V3-2iV| 3iVc-27V| 7Vc(2Ar2_3) 3 - 2iV2 

, „ , 2(iV,^ + 1) , 2(jV,^ + l) X 

The element Vi is identical to the tree-level vertex, whereas V2 and Vz represent the two only additional Bose 
symmetric structures that can be built from eq.®. The vertex is then represented by 

''Kti ipi,P2,P3) = E ''r. (Pi,P2,P3) (11) 

The object F4f(p2) appearing in the running coupling ^ is then related to (|lip by 
p4g^p2-j _4g (pi,P2,P3), whcrc all external scales ^ P2 P3 ^ Pi ■ P2 ^ Pi ■ Ps ^ P2 ■ P3 ^ P^ ■ We will 
come back to this coupling in section fVl 
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B. The DSE for the four-gluon vertex 



Having constructed a suitable representation of the four-gluon vertex we now discuss the structure of its Dyson- 
Schwinger equation. In compact notation this equation reads [soj : 




where all internal propagators are to be understood as fully dressed and the shaded circles denote reducible vertex- 
functions. The decomposition of these functions into proper (i.e. one-particle irreducible) vertices is given in appendix 
[Cl Here also higher n-point functions (n = 5, 6) appear which satisfy their own Dyson-Schwinger equations. Since 
in general one cannot solve the resulting infinite tower of equations at once, we have to truncate the vertex-DSE, 
eq. (jl2p . in a physical reasonable way. The truncation scheme that will be applied is as follows: 

• The four-gluon vertex will be reduced to the subset of structures discussed above. In particular we use the 
Bose-symmetric representation, eq. (jlip . 

• The fully dressed ghost and gluon propagators in the internal loops are taken from their own coupled system 
of DSEs. These have been solved in 1221 without taking into account any effects of the four-gluon vertex. By 
comparison with lattice calculations (32), 'ssj one finds that this approximation in the propagator DSE leads to 
errors of the order of ten percent in the mid- momentum region only 20]. The far infrared and the ultraviolet 
are unaffected. We therefore employ these solutions in this work and leave an inclusion of the back-reaction of 
the four-gluon vertex on the propagators for future studies. 

• Due to the complexity of the four-gluon vertex DSE it seems justified to reduce the number of diagrams contained 
in our investigation to the ones that give dominant contributions in the infrared and ultraviolet momentum 
region. Since these limits are under analytical control (see section Hill and refs. we can identify these 
diagrams safely. In the infrared, the leading diagram is the ghost-loop (b), whereas in the ultraviolet leading 
contributions can be expected from all one- loop diagrams. 

• These diagrams then contain higher n-point vertices, that will be reduced to two- and three-point functions 
using a skeleton expansion (i.e. an expansion in full vertices and propagators). 

• Selfconsistency effects of the four-gluon vertex will be neglected, i.e. we drop all diagrams on the right hand side 
that contain the four-gluon vertex (e.g. the diagrams (c),(d),(e),(f) in eq. (H^. While this approximation greatly 
reduces the complexity involved in the numerical treatment of the DSE it does not affect the infrared behavior 
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of the resulting four-gluon vertex, since the ghost loop (b) is the dominant diagram for small momenta (c.f. 
above). In the ultraviolet momentum region, however, this omission leads to a one- loop running of the vertex 
not in agreement with perturbation theory. We remedy this drawback by the use of an effective three-gluon 
vertex in diagram (b). 

• This effective three-gluon vertex obeys the correct IR power-law and generates the correct UV behavior of the 
four-gluon vertex under absence of the diagrams (c),(d) and (e). A similar effective construction has been used 
previously in the DSEs for the ghost and gluon propagators, where results close to corresponding ones from 
lattice calculations have been obtained [2^ . 

• The dressed ghost-gluon vertex will be replaced by the bare vertex. This approximation is well justified not 
only in the ultraviolet but also in the infrared momentum region. This property has already been conjectured 
by Taylor in the early seventies [3^ and has recently been verified numerically in continuum as well as lattice 
calculations 'ss', 's^, [s^] . It also agrees with the all-oder analytical analysis of the DSEs performed in [l^, [3 , 
cf. section UlI Bl 

The resulting approximation of the Dyson- Schwinger equation of the four-gluon vertex then reads 




(13) 

Here 'perm.' denotes permutations of the three external dressed legs of the ghost-box and the gluon-box diagram. 
The subscript 'symm' indicates that we average over all possible locations of the bare vertex in the diagrams thus 
restoring Bose symmetry by hand. The dressed ghost-gluon vertices are taken bare and the dressed three-gluon vertex 
is given be the following ansatz: 

Z((72)( — ) Z(p2)(— ) ^ > 

with the one-loop anomalous dimension of the ghost S ~ — ^ and the dressings G{p^) and Z{p^) of the ghost and 
gluon propagators. The symbol ^p^F denotes the bare three-gluon vertex. This ansatz preserves the correct UV 
anomalous dimension of the full four-gluon vertex, as well as the correct IR power law of the three-gluon vertex in 
the scaling scenario reviewed in the next section. 



III. ANALYTICAL RESULTS IN THE INFRARED AND ULTRAVIOLET MOMENTUM REGION 



A. Multiplicative Renormalizability 

The truncation of the DSE for the four-gluon vertex, eq. (jlSp . is given explicitly by 

^ abed — ^4 (o)J- abed 

\ J \ ^7 / \perm.]symm 

+ (z^g^ j J^^SYP'PP" Dl<I,^,T^'^^" Dll^,T^^'^^^^" D^,.,, l^^T^'"''" p,Tabc^ (15) 

\ J \ / \perm.;symm 

where the color factors have been subsumed in a factor Tabcd — ifb'aa' f&bb' fd'cc' fa'dd') (recall that we assume the 
ghost-gluon and three-gluon vertices to be proportional to /"^"^ and the propagators to be diagonal in color space) . The 
symbols (o)!"' and (q^T"^ and ^q^T^^^J denote the bare ghost-gluon, three-gluon and four-gluon vertices respectively. 
All momentum arguments have been omitted for brevity. Note that this equation is already divided by a factor g"^ 
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coming from the full vertex on the l.h.s. of eq. (|13p. Before we embark in the analytical analysis of this momentum 
dependence, we wish to show that this truncation scheme preserves multiplicative renormalizability of the four-gluon 
vertex DSE. 

To this end we need the relations between the renormalized and the unrenormalized but regularized Green's functions 
of the theory. The former ones are functions of the renormalization point /i^ (in addition to their momentum 
dependence), whereas the latter ones depend on the regularization scale. If the regularization is performed by a 
momentum cutoff A these relations are given by 

gi^i)Zg{ii,A) = .g(A) (16) 

DGip,ti)Z3i^x,A) = Dg{p,A) (17) 

Dp^{p,fi)Z3ifi,A) = Dp,ip,A) (18) 

^OTip.^fiY^^^'' = Zi{fi,A)^sr[p^,A)P^^^'' (19) 

^^Tip,, ^i)P^^ = Z,{p,Afsr{p^,A)P^^ (20) 

Tip,,^^r = Z,{ti,AfaTip,,Ar (21) 
They are complemented by the Slavnov- Taylor identities 

Zi — ZgZ^^^, Zi — ZgZj,z\J^ ^ ZiF = ZgZ\l''Z2^ ^4 = ^3 ' (22) 

One can then analyze the dependence of the ghost and gluon box diagrams on the renormalization point /i of the 
theory. We obtain for the renormalization point dependence of the ghost box diagram g}J^^^cd^^^) 

.'.rl7(M') - |^|()^[^i(m')]'||-^ = {Zg{i?)\\Z,{i?t = , (23) 

and for the gluon box diagram giX'^^^d^^) 

uKtl^l^^) - - {Zg{i?)\\Z:,{^?)f = Z,{i?) , (24) 

As a result, all diagrams are proportional to Z4(/i^), which guarantees the multiplicative renormalizability of the 
vertex DSE in our truncation scheme. 



B. Yang-Mills Green's functions in the infrared 



The infrared behavior of the four-gluon vertex can be determined from its Dyson-Schwinger equation by means of 
analytical techniques. Before we demonstrate the details of such an analysis we have to shortly summarize previous 
results on the infrared scaling of general one-particle irreducible Green's functions of Yang-Mills theory. 

The basic idea, followed in [1^, to determine the infrared behavior of one-particle- irreducible (IPI) Green's functions 
is to investigate their Dyson-Schwinger equations order by order in a skeleton expansion {i.e. a loop expansion using 
full propagators and vertices). The analysis rests upon a separation of scales, which takes place in the deep infrared 
momentum region. Provided there is only one external momentum scale p^ « Aqcd much smaller than Aqcd, a 
self-consistent infrared asymptotic solution of the whole tower of Dyson-Schwinger equations for these functions is 
given by 

r"'"(/) - (p2)(«-™)K_ (25) 

Here r"''"(p2) denotes the dressing function of the infrared leading tensor structure of the IPI-Green's function with 
2n external ghost legs and m external gluon legs. The exponent k is known to be positive [13, H^. 

A special instance of the solution ((25|) are the inverse ghost and gluon dressing functions T'^--^{p^) = G^^{p'^) and 
p0,2(p2-j _ z~-^{p'^), which are related to the ghost and gluon propagators via 

D-ip^)^-^, D,Ap^)=(s,.-'-^)^. (26) 

The corresponding power laws in the infrared are 

G{p')^{pT\ Z(p2) ^ (p2)2. _ (27) 
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For a bare ghost-gluon vertex in the infrared, justified by lattice calculations (Sy, |37| and also in the DSE-approach 
[sst . one obtains k = (93 — -\/1201)/98 « 0.595 [13, Ull- Possible corrections by regular dressings of the vertex in the 
infrared have been investigated in [l^, where an interval 0.5 < k < 0.7 has been given. Thus, although the precise 
value of K is hitherto unknown and depends on the truncation scheme, the variation is quite small and not important 
for the results presented in this work. 

An interesting consequence of the solution (|25p is the qualitative universality of the running coupling in the infrared. 
Renormalization group invariant couplings can be defined from either of the primitively divergent vertices of Yang- 
Mills-theory, i.e. from the ghost-gluon vertex [gh — gl), the three-gluon vertex (3g) or the four-gluon vertex (4(7) 
via 

a3'^-9\p^) = ^G2(p2)z(p2) const/TV,, (28) 

47r 

47r 

47r 

Using the DSE-solution ([25| it is easy to see that all three couplings approach a fixed point in the infrared. This 
fixed point can be explicitly calculated for the coupling (j30p . Employing a bare ghost-gluon vertex one obtains 
a9h-gi{Q) « 8.92/iVc [3]. 

We emphasize that the eq. (fSS]) solves the untruncated system of DSEs and the corresponding equations from the 
functional renormalization group. Thus, although k depends on a truncation scheme, ([25]) does not. It is furthermore 
the only possible solution of both systems in terms of irrational power laws T?] . The resulting behaviour of the gluon 
and ghost propag ators agree well with the predictions deduced in the Gribov-Zwanziger and Kugo-Ojima confinement 
scenarios [23, [33]. Nevertheless there is a caveat here: lattice Monte-Carlo simulations have not yet been able to 
verify the relations ([?f|) . In fact, very recent results on large lattice indicate that the exponent of the gluon dressing 
function ma y b e close to k « 0.5, whereas the corresponding value for the ghost dressing function may be considerably 
smaller [40|, l4i|. These findings allow for at least two possible interpretations: they may indicate a different infinite 
volume limit than expressed by eq. ([25|) . or they may be attributed to Gribov-copy effects associated with gauge fixing 
on large lattices. General considerations on the confining properties of QCD suggest the latter interpretation 
Pending further clarification we will therefore employ the behavior eq. (|25p for the purpose of this work. 



C. Infrared analysis 

1. The ghost box diagram 

According to the general analysis of , the infrared behavior of the four-gluon vertex in the presence of only one 
external scale is given by 

r°-4(p2) ^ (p2)-^-, (31) 

see eq. ()25p above. This solution is generated by the ghost contributions to the vertex-DSE, i.e. in our truncation 
by the ghost box diagram. In the following we will verify this result for one particular momentum configuration and 
determine the corresponding coefficient of the power law. This will be useful for two reasons. First the result provides 
a welcome consistency check to our numerical calculations. Second, and much more important, together with the 
corresponding result for the gluon propagator it will give us the value for the infrared fixed point of the running 
coupling from the four-gluon vertex. In principle, the running coupling can be calculated for every basis component 
projection of the full four-gluon vertex. However, matching with perturbation theory in the ultraviolet momentum 
region demands to perform this analysis with the tree-level tensor structure, which will be done in the following. 
The particular momentum configuration we choose for our analysis is given by 
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It has the merit that it is invariant under permutations of the three dressed legs. Thus all permutations give the same 
results which can be taken into account by a factor of six in front of the integral. ^ 

With bare ghost-gluon vertices and projected onto the tree-level tensor the ghost box is then given by 

, 2. g^N, f , 2 2,^, G{q + p) G(q + 2p) G{q + 3p) G{q) 

OO TT 

- I Idq^dep^qU^-ie f^^^^^'^^^^^^^''^^^'^^''^^^ 



36(2^)3 J J ^ ^ ^ ' ' {q+pY {q + 2pY {q + Sp^ q" 





(32) 



Note that this contribution is already Bose- symmetric in our truncation scheme with bare ghost-gluon vertices. The 
factor six from permutations of the external legs is already included here. 
Since the internal ghost dressing functions are infrared divergent, i.e. 

G{p)=B{p'y\ (33) 

for p^ « Agpjj the integral is dominated by loop momenta where the internal momentum is of the same order as 
the external scale p^ . We can thus replace the internal ghost dressing functions by the infrared asymptotic expression 
This leads to 

OO TT 

^«r,, ^ j J dqHOp^q' sm\6){q + p)-2(i+'^)(q + 2p)-2(i+-)(q + ip)-^i^+-) q-^i^+-) . (34) 



We then divide the momentum integration range into three parts from [0,^^], [p^ , IQp^] and [lOp^, oo] and denote the 
corresponding contributions by /a,/b and Ic- 

Of course, replacing the internal ghost by the infrared asymptotic expression (j33p would be a poor approximation 
if the contribution Ic were to dominate the total integral I — la + h + Ic- However, this is not the case. Ic can be 
evaluated using a Taylor expansion and we find its contribution to be extremely small compared to la + lb provided 
the lower bound of this integral is chosen large enough. This is indeed the case for our choice q^ > lOp^ and we may 
therefore neglect Ic- 

To evaluate the first integral, la, the approximation 

(g + p)'('? + 2p)2(g + 3p)2«36/ fl+ +2^cos(0)^ (35) 

\ \apj ap J 

is employed, with a parameter a > 1. This parameter can be determined numerically; we find a ~ 1.886 and obtain 

^ g^NcB^ ) 2 P^q' / ,„ sin>^ 

"~ 36(27r)3 7 '^'^ (36pV)"+W / / ^2 x3(.+i)- 

- +2^cos(^) 



The angular integral can be evaluated with eq. (|Dip . yielding 



g^NcB^ } ^ ^ p^q'^ ^ f ^ l\ ^ /^o... , -, ^ o.. , i . o. 9 ^ 



" -§6(2^7'^'^ (3##TT^^2'2j^^ I- 



^ At first sight one may believe that an even simpler kinematical choice is possible, namely pi = P2 = — P3 = — P4- However, we found 
that such configurations lead to results which are not stable with variation of an numerical infrared cutoff e. This indicates, that the 
emergence of such a kinematic situation as a limit of a more general setup is not free of singularities. Such 'soft' or 'coUinear' singularities 
arise in addition to the 'overall' singularity 1)31^ of the four-gluon vertex. In this work we will not touch upon these soft singularities 
and leave this issue for future studies. 
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Abbreviating z ~ ^ one then obtains ^ 



1927r (36)'-+i ^ 



a{n'^)NcB^ 1 r(2-K) 
192^? (36)«+i r(3 - k) 



( 3(k + 1), 3k + 1; 3; — z 



Fa 2 - K, 3(k + 1), K + 1; 3 - K, 3; ^ x (p^) 



(38) 



with a^M^) — g"^ /{Att). The last integral has been solved with the help of eq. (|D2p . Inserting n = (93 — vT20T)/98 
(cp. the text below eq. (HZ])) and a — 1.886 one finds 



la « 9.49 • 10-^ • a(M') • TV, • • (p^) 



2\-4k 



(39) 



which agrees with the power counting analysis, eq. pip . 

Now only the part /f, where the loop momentum is of the same order of magnitude as the external momenta is left. 
It can be evaluated using a Chebyshev-expansion (in the loop- momentum and the polar-angle) , see appendix ID 21 for 
details. Renaming variables as x — , y = and abbreviating 



x{9x + y + Gy^Ey cos{e)y'-^^''^ 
g{y,9) - sin4((?)(l + j/ + 2Vycos(0))-^'+"^(4 + 2; + 4V^cos(0))-(i+^) 
X (9 + y + 6^cos{e))-^^+''^y-^^+''\ 



and 



-I cos 



(fc - l/2)7r 



N 



1 



yi = xyk^ X \ - cos 



(/-l/2)^ \ 11 
N' 2 



one finds 



(40) 
(41) 



(42) 



'' 36(27r)3 N N'^^ ' 



N' 



El / /, \ cos(z7r) + 1 



fc=i \ i=i 

iV-l 



E cos(j7r) + 1 



1=2 



cos 



i=2 

N' 



k=l 
N'-l 



1=1 

i(fc- 1/2)77 



i{l - l/2)7r 



, /, X cos(i7r) + 1 

E5(y^^fc)+ E i_-2 E 



N' 



. 1=1 



i=2 



1=1 



l/2)7r 

W' 



(43) 

This expression can be evaluated numerically. It turns out that the expansion is well converged with N = N' = 20. 
This yields 

h « 9.49 • 10"^ • a{^i^) -Nc-B^- {p^y^". (44) 
This contribution is almost exactly a factor of ten larger than la- 



Actually eq. I I38I I directly shows the abovementioned soft singularity occurring when the momentum configuration pi = P2 = — P3 = —Pi 
is chosen. In this case the approximation eq. I|35| l becomes exact and a = 1. However the hypergeometric function in eq. I|38| l does only 
converge when \a\ < 1. 
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Putting all pieces together one finally finds 

'"^TgH [p") = la+h+Ic 

« 1.04 • 10-4 AT, .B4.(p2)"4«_ (45) 

This result will be used in section |Vl where we discuss the infrared behavior of the running coupling. 

2. The gluon box diagram 

The IR-behavior of the gluon box diagram can be estimated by power-counting using eq. (j25p . In this diagram, 
there are four gluon propagators along with three three-gluon-vertices. The three gluon vertices behave like 

^«^3,(p2) = c•(p2)-'^ (46) 

when 0. Together with the four gluon propagators, which contribute to the IR divergence like C ■ (p^)^'" one 

gets for the gluon box diagram 

^«^,, = c"•(p2)-^ (47) 

The gluon box thus is only subleading in the infrared in agreement with our general considerations in section [III Thus 
the coefficient C" is of only minor interest and will not be computed here. The power law behaviour (p^)"" is well 
reproduced by our numerical results for the gluon box. 



D. Ultraviolet analysis 

1. The ghost box diagram 

It is known from resummed perturbation theory that in the ultraviolet momentum region the dressing function of 
the ghost propagator can be described by the asymptotic expression 

Gip')^Gi^i')(^cJ\og(^^^+l^ , (48) 

with the one-loop anomalous dimension 6 — —9/44, tu = llNca{fj,'^) /12Tr and some renormalization point /i^. This 
behavior is reproduced by the Dyson-Schwinger equations for the ghost propagator Q. Plugging this into eq. ([32]) 
and using dimensional regularization we arrive at 

r.* (P ) - - 36,2,)J J ''"PI »n («) ip7^2 f^^^2 



;.,o,(i^).l) (.■os(^).l^ 



5 



which describes the ultraviolet behavior of the tree-level projection of the ghost box dressing function. Since the 
diagram is dominated by the region where the loop-momentum is larger than the external momenta, it is justified 
to employ an angular approximation: all arguments of the logarithms are replaced by the loop momentum q^. The 
integration interval can then be restricted to [p^, oo]. Furthermore, the denominators are approximated using eq. (j35p . 
After evaluating the two trivial angular integrals and the third one using ea. (jDl[l we obtain 



^^''(^ ) - - 2592 (2 V FMf) > i ) I'' ' ' 2^ 2 + ^' ^ j l^^^^ If J + ' 

(50) 

The hypergeometric function has the series representation 

2Fl(a,/3,7,^)=f:^^^^^ (51) 
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with the Pochhammer symbol {a)j as introduced in appendix ID3I The remaining integral can be evaluated with the 
help of eq. (|DT6|) . Wc find 

where the dj are different regulator dimensions, one for every order j of the expansion eq. (|5ip . The divergence of the 
integral is absorbed into the coefficients of the logarithm and finally cancelled by the renormalization procedure. In 
order to match our results from the numerical calculations, where a momentum cut-off regularization will be employed, 
the renormalization condition for the analytical result is chosen such that the numerical and the analytical results 
agree at the renormalization point /i^. 

As can be seen from the Slavnov- Taylor identity Z4 = Z^/Z^ the anomalous dimension of the four-gluon vertex 
in the ultraviolet momentum region should equal —7 + 26 — 1 + 46 = 8/44, where 7 = —13/22 and 6 = —9/44 are 
the anomalous dimensions of the gluon and the ghost propagator and 1 + 2(5 + 7 = 0. The result 46 ~ —36/44 found 
here is negative and leads to a vanishing contribution in the ultraviolet. Thus the ghost box is subleading at large 
momenta and the leading contributions have to come from the gluonic diagrams. 



2. The gluon box diagram 

With the abbreviations Qq = q^, qf = [q f?! = + and q^ = {q + 3p)^, the gluon box integral reads 

f .^Z{ql)Z{ql)Z{ql)Z{ql) 



-r^KXau ( 2\ 9 f d q , 2 2 , 

abed [P ) = J J^^yP '9 



2 9 9 9 

Ti ii % % 



(Z(gg))(^) \{Z{q^))^^^ ) \ {Z (qi))'^^^ 

with the model three-gluon vertex given in eq. (I14p . The kinematic kernel K (p^, g^, 0) stems from Lorentz contractions 
after projection on the tree-level vertex. This kernel is complicated and lengthy and we therefore omit its explicit 
form. It has the general structure 

6 10 

K{p',q',e) = ^cos"0 X J2 (P'f (9')'"^ ■ (54) 

n—Q m— 

— 

For even n we can replace cos'^ 9 — (l — sin 6*) ^ . If n is odd, one has to factor out one cosine. Then the kernel has 
the structure 



,„ ^ ™ f (1 - sin^ 6*) neven 
,i=Om=o I cost^(^l — sm 0) nodd 



This form allows using eq. (|D2[) and (|D3p for the analytic calculations. Similar to the ghost dressing function the 
ultraviolet behavior of the gluon dressing can be written as 

Z(p2) = Z(/i2)(^^log(^^)+iy. (56) 

As before in the ghost box we employ the angular approximation qo — qf — q2 ~ Q3 — T' ij^i the logarithm, use the 
approximation eq. (|35p and dimensional regularization. We then find 



1+45 



r.^ [p' 



/\uj log -T- -^ i 
dUK{p^,q\e) ^ \!l2_J. ^. (57) 

36P«9^ (1 + (i) + 2^ cos(0) j 
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Evaluating the trivial angular integrations and again restricting the integral on the interval [p^ , oo] then leads to 

2^(^)r(^) {G{f.')y'-'' 



^^'^^^ 72(2^)'^ T^^) {Z{t,^))'+'' 

^ 1+4(5 

X [dy I d9 i,^)^sin^'^e ^ ""^'^^'^'^ f.logf^) +l) (58) 

{ 36.3(l + ^ + 2V^cos(.)) ^ ^ 

with eqs. (|D2|l . (|D3|l and (|DT6)) one then finds 

/ / 2 \ \ ^^^^ 
^^r,Kp') = C • (^^log(^^j +lj , (59) 

with a regularized factor C. This divergent factor is extremely lengthy and we therefore refrain from giving it 
explicitly here. However we note that the divergence is such that it matches the one of Z4 in the four-giuon DSE thus 
guaranteeing a finite result for the four-gluon vertex on the left hand side of the DSE. The momentum dependence of 
(|59p and in particular the anomalous dimension 1 + 45 of the logarithm is in agreement with the expectations from 
the Slavnov- Taylor identity for the four gluon vertex renormalization factor as discussed above. 



IV. NUMERICAL RESULTS FOR THE FOUR-GLUON VERTEX 



A. Numerical methods 



For the numerical investigation the subtracted version of eq. (|13p . is considered. 




with some subtraction point which is set equal to the ultraviolet momentum cutoff = 10^" GeV^ for reasons 
of numerical stability. Furthermore we also introduce an infrared cutoff e^, which is chosen as — 10^^" GcV^. The 
numerical integration is carried out using a Gaufi-Legendre algorithm on a logarithmic grid. 

In eq. ([50]) . the renormalization constant Z4 stemming form the tree-level vertex after projection, is replaced by an 
integral and the full dressed vertex function at the subtraction point. Its value has to be fixed by renormalization. 
This will be done such that the running coupling from the four-gluon vertex has the same value as the running coupling 
from the ghost-gluon vertex at the renormalization point /i^. As input we use the value a{fi^ = 1.713 GeV^) = 0.97 
determined in [22| within a momentum subtraction scheme. 

Before we present our results we need to discuss two further technical points: 

• We explicitly checked that our results are independent of the ultraviolet and infrared cutoffs. This is indeed the 
case, if these cutoffs are at least three orders of magnitude larger/lower than the largest/lowest of the external 
momenta. In addition, due to the complicated kinematics of the gluon box we had to use quite a large number 
of sampling points for the radial momentum integral (typically 5000 points on a logarithmic grid) . To further 
improve the numerical accuracy the momentum integral has been split into three parts, integrating the infrared 
up to a small region ± Ap^ around the external scale, the small region itself and then up to infinity. In 
the numerical calculations Ap^ is chosen as Ap^ = O.Olp^. In the numerical angular integral we find that is is 
numerically advantageous to integrate over the cosine of the angle. 

• The crucial assumption in our infrared analysis of the vertex-DSE was, that the integrand of the diagrams are 
dominated by loop momenta of the order of the (small) external momentum. This we verified explicitly also 
numerically. 
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FIG. 1: Results for the full four-gluon vertex (top) and the ghost box diagram (bottom left) and the gluon box diagram (bottom 
right) in the kinematical section pi — — 3p,p2 = ps = P4 = p. The results for the ghost and gluon boxes include the symmetry 
factors. The three curves in each diagram correspond to projections onto the three Bose-symmetric structures given in eo. 1101 



As input for the ghost- and gluon propagators we take the following analytical expressions 

a(x) = ■ R(x) = ^ ' ^ ' (61) 

He + a^{x/^l^^r^^h,[x/Kl^^f^y I + c[x / Kl^^Y + d{x/ Kl^^Y^ ^ ^ 



Z(x) = 



a{x) 



1+25 



R\x)- 



G{x) = 



a{x) 



R-\x), 



(62) 



with parameters 



a(0) 




ai 


02 


bi 


62 


c 


d 




8.915/iVc 


0.97 


1.106 


2.324 


0.004 


3.169 


1.269 


2.105 


0.714GeV 



and the anomalous dimension S = —9/44 of the ghost. These expressions have been fitted to the numerical results of 
[2^ 1 for the coupled system of DSEs for the ghost- and gluon propagators. 



B. Results 



We first present numerical results for the specific kinematical situation given by 
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which matches the one used in our infrared and ultraviolet analysis. The results are presented in Fig. [T] On the 
top diagram we display the full four-giuon vertex in this kinematical setup projected onto the Bose-symmetric tensor 
structures given in eq. p^ . Recall that the structure Vi is identical to the one of the bare four-gluon vertex ('tree-level 
structure'). The results match nicely our expectations from the analytical analysis in section Fill CI All structures 
of the four-gluon vertex diverge like (p^)"**" in the infrared. This divergence is driven from the ghost-loop diagram, 
as can be seen from comparing the full result with the contributions from the ghost box and gluon box diagrams 
displayed in the lower panel of Fig. [1] Concerning the infrared coefficients it is not the Bose-symmetric tree-level 
structure that dominates but one of the non-tree-level counterparts. 

All curves show a characteristic scale of a few hundred MeV, where the infrared power law behavior bends towards 
the logarithmic, perturbative behavior in the ultraviolet momentum region. Certainly, the magnitude of this scale is 
inherited from the input ([62|l for the ghost and gluon dressing functions and represents the scale Ay m of Yang-Mills 
theory generated by dimensional transmutation. 

In the ultraviolet momentum regime we also reproduce the analytic behavior of the tree-level structure determined 
in section IIII CI Here the leading contribution stems from the gluon box diagram. Similar to the infrared, we also 
observe that the non-tree-level structure V2 has the largest coefficient of the three structures considered. However, 
this will change for even larger momenta, since the logarithms appearing in V2 and V3 have negative anomalous 
dimensions, while the tree-level structure Vi has the correct and positive anomalous dimension 1 + 45 in agreement 
with resummed perturbation theory. 

At first sight, it seems counter-intuitive that the structure V2 dominates the vertex also for the relatively large 
momenta considered in our calculations. However, this dominance has a natural interpretation: it is the prefactors 
of these contributions stemming from the corresponding color contractions that give large relative coefficients of the 
order of 10^ between the Vi and V2 projections. These also appear in first order perturbation theory, i.e. with no 
internal dressings from propagators and vertices in the ghost and gluon box diagrams. We explicitly checked that the 
relative ordering of the contributions Vi,V2 and V3 is the same in this case. This shows that the ordering appearing 
in Fig. [1] is not an artefact of the truncation of the ghost-gluon and thrce-gluon vertices. The relative magnitudes of 
Vi, Vi and V3 may however be modified by the inclusion of the missing one- loop diagrams (c),(d) and (e) of Eq. (|12p . 

In Fig. [2] we also present a calculation for a different kinematical situation with the two independent Lorentz 
invariants and P2 ^-iid p\- P2 — |pi||p2|, Pz — Pi, Pi — ^Pi — 2p2 (all four momenta P1..4 are defined to flow into 
the diagram), i.e. 



As can be seen from Fig. [2] we find an infrared divergency when all momenta go to zero with the power law (p^)"^'' 
satisfied in the presence of only one external scale, i.e. in a cone around the diagonal of the diagram. The behaviour 
of the vertex dressing function for kinematics at the edges of the diagram is nontrivial and may indicate additional, 
weaker kinematical singularities present when one or more external momenta are held fixed. Note that the numerical 
problems in these case are quite intricate, since the presence of an additional scale involves huge cancellations for 
some kinematical points. Here we dealt with these problems by employing an adaptive framework for the angular 
integration routines. Since the main focus of this study is the running coupling defined along the diagonal of Fig. [2] 
we postpone further discussion of the general kinematical behaviour of the vertex to future work. 



P2 




(-Pi - 2^2) 



15 




FIG. 2: Results for the ghostbox in the kinematical situation with pi • P2 ~ |pi||p2|, P3 = P2 and p4 = — pi — 2p2 (aU four 
momenta pi.A are defined to flow into the diagram). 



V. THE RUNNING COUPLING 



The running coupling from the four-giuon vertex has already been given in eq.Q and is repeated here for the 
convenience of the reader: 

Here r*3(p^) denotes the dressing of the tree-level vertex structure Vi and Z{p'^) denotes the dressing function of the 

2 

giuon propagator. The value — ) = 0.97 has also been given before. The resulting momentum dependence 
of the coupling is shown in Fig. [3] along with the running coupling from the ghost-gluon vertex. In the ultraviolet 
momentum region we observe that both couplings run like the usual inverse logarithm as well known from perturbation 
theory. Here we have universal behavior as dictated from gauge invariance. In the mid-momentum region we observe 
a steep rise of both couplings up to values of order a ~ 1. The coupling from the ghost-gluon vertex then keeps rising 
before freezing to an infrared fixed point, whereas the coupling from the four-gluon vertex decreases dramatically 
until it reaches a very small but nonzero fixed point in the deep infrared. 

Before we discuss the implications of this behavior further, we wish to verify this numerical result from our analytical 
calculations in section UlI CI There wc found that 

'"rg^{p') « 1.045- 10-4 •a(Ai2).Ar^.B4.(p2)-4. 

= ■ ain^) ■ ■ ■ {py'\ (64) 
In the infrared the gluon dressing function obeys the power-law 

Z{p^)^A{p^f\ (65) 

We then obtain 

a,, -> 0) = a (m^) • C,, ■ aip') -N^-B^- {p^y'^ ■ [a [p^f 

= [a{p'')-A-B^f-N,-Cgh, (66) 

This expression is manifestly RG-invariant, since we know from the coupling of the ghost-gluon vertex that aU^AB"^ 
is an RG-invariant. Also, since q;(/z^) = (7^/(47r) ^ l/Nc and A and B are separately independent of Nc [H, [l^l as is 
Cgh, the coupling is proportional to 1/Nc in agreement with the large Nc counting rules. 
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FIG. 3: The running coupling from the four-gluon vertex compared to the coupling from the ghost-gluon vertex from ref. [22 



Note that the combination a{^^)AB^ is equivalent to the running coupling from the ghost-gluon vertex at zero 
momentum [l3|. Thus we can rewrite eq. ([M)) as 



With aS^-s^iO) « 8.92/iVe [M] and Cgh = 1.045 • IQ-* we then obtain 

0.0083 



a4g {p O) 



N, 



(67) 



(68) 



This value agrees well with our numerical result. 

In order to assess the implications of our findings it is worth noting that the corresponding coupling from the 
three-gluon vertex has a very similar behavior than the one of the four-gluon vertex shown in figure [3l In particular 
it also has an extremely small fixed point in the infrared, a maximum at intermediate momenta and a perturbative 
logarithmic tail in the ultraviolet (4^ . We therefore find universality in the ultraviolet momentum region, as required 
from gauge invariance. In the infrared, all three couplings are qualitatively similar in the sense that they all go to 
an infrared fixed point (as already emphasized in [15]). However, there are huge qualitative differences between the 
coupling involving ghosts, agh-gi, and the other two couplings, asg and a^g, that only involve gluonic correlators. 

In this respect it is important to note, that the smallness of the infrared fixed point of the four-gluon vertex is 
rooted in the structure of the vertex-DSE. In section (jlll C ip we found, that the diagram involving ghosts is the one 
that gives the leading infrared behavior and determines the coefficient of the infrared power law Cgh{p^)~^'^ of the 
four-gluon vertex. One reason why the coefficient Cgh is small is a factor 1/216 stemming from the projection on the 
tree-level tensor structure. The four propagators in the loop generate further suppression. Therefore the smallness 
of this coefficient can be attributed entirely to the structure of the DSE and does not depend on our choice for the 
ghost-gluon and three-gluon vertices. 

Of course, the ghost box is only the first term in the skeleton expansion of the original ghost related diagram in 
the full DSE, eq. (jl2p . It is known from ref. 15] that all terms in this expansion share the same infrared power law 
behavior and will therefore contribute to the coefficient Cgh ■ Certainly, we cannot exclude that a summation of these 
terms will result in large changes compared to our value of Cgh- 

However, there is a systematic argument that indicates this may indeed be not the case. Consider once again 
the Dyson-Schwinger equation for the gluon propagator. In Fig. [1] we compare this equation with the corresponding 
flow equation in the functional renormalization group framework. From a systematic point of view, the DSE can be 
viewed as an integrated flow equation, so the physical content of the two equations is the same. However, there is 
an important structural difference between the two equations: In every diagram of the DSE we have one bare vertex, 
whereas in the flow equation, all vertices are dressed. As has been noted in this leads to an interesting situation 
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FIG. 4: The Dyson-Schwinger equation (top) and the functional flow equation (bottom) for the gluon propagator. In the flow 
equation, crosses denote insertions of the infrared regulator, which cuts off the theory at or around a scale k. 

for small momenta: In the DSE the ghost-loop (a) is the only diagram responsible for the leading contributions in 
the infrared, all other diagrams are suppressed by powers of momenta. In the flow equation, however, all diagrams 
share the same infrared exponent and therefore contribute to the coefficient of the resulting power law for the gluon 
dressing functions. Now, both equations are exact as they stand, so they should give the same results in the infrared. 
We also know that the ghost-gluon vertex is almost bare in the small momentum region. For this reasons we expect 
that the ghost-loop diagrams in the flow equations (first line) should be roughly similar in effect to the ghost-loop 
(a) in the gluon-DSE 43,, 44]. Thus the infrared coefficients of all other diagrams in the flow equations have to either 
cancel, or should be much smaller than the coefficient of the ghost diagrams. Our result for the coupling from the 
four-gluon vertex together with the tentative result for the three-gluon vertex [43| indicates exactly this: since all 
these diagrams are roughly proportional to their corresponding coupling the gluonic diagrams in the flow equation 
are parametrically suppressed compared to the ghost-diagrams due to the smallness of the gluonic couplings in the 
infrared. This offers a natural explanation how the DSE and the renormalization group framework can agree in the 
infrared. In turn, this point suggests that the smallness of the four-gluon coupling may indeed be an effect which is 
robust beyond the leading order in the skeleton expansion. 

In fact there is a further argument supporting this scenario. In [23 | Zwanziger gave good arguments for an infrared 
effective theory dominated by the Faddeev-Popov determinant. He argued that all purely gluonic interactions switch 
off in the infrared and it is the geometry of the gauge group which then controls the infrared dynamics via the ghost 
content of the theory. This is exactly what we found here. 



18 



VI. SUMMARY 

In this work we investigated the nonperturbative structure of the four-gluon vertex from (a truncated version of) its 
Dyson-Schwinger equation. We identified analyticaUy the leading infrared and ultraviolet terms of this equation and 
found good agreement of this analysis with our numerical solution. We investigated the behavior of the three Bose 
symmetric tensor structures that can be constructed from a subset of the complete tensor basis of the vertex. The 
dressing functions of these three structures all show an infrared singular behavior with power laws in agreement with 
the results from naive power counting [Tsj . In the ultraviolet momentum region our solutions reproduce resummed 
perturbation theory. 

The central result of our work concerns the running coupling from the four-gluon vertex, built from a combination 
of vertex dressing and the dressing function of the gluon propagator. Although in the ultraviolet momentum region 
the coupling agrees nicely with the one from the ghost-gluon vertex (as it should, according to gauge invariance), in 
the infrared we observe strong deviations. Whereas the coupling from the ghost-gluon vertex develops an infrared 
fixed point at around agh-gi{0) ~ ^/Nc, we find a much smaller fixed point at around a3g(0) « 9 • 10~^/Nc for the 
coupling from the four-gluon vertex. 

Certainly, the stability of this finding has to be checked wrt further improvements of our truncation scheme. These 
have to include a study of the two-loop diagrams in the skeleton expansion of the ghost-part of the vertex DSE, since 
(only) these terms have the potential to change the infrared coefficients of the vertex and therefore the value of the 
infrared fixed point. However, on general grounds we are confident that the smallness of the running coupling from the 
four-gluon vertex is an important property which is stable wrt these improvements. As discussed in the last section, 
the reason is that this fact explains why the Dyson-Schwinger and the functional renormalization group equations 
for the two point functions of Yang-Mills theory agree in the infrared, although their structure is quite different. 
With small couplings from the three- and four-gluon vertices gluonic contributions to the infrared behaviour of the 
ghost and gluon FRGs are parametrically suppressed, leading to ghost dominance in agreement with the results from 
the DSEs. This finding also supports the notion of an infrared effective theory dominated from the Faddeev-Popov 
determinant proposed in (24| . 
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APPENDIX A: THE TENSOR BASIS 

The tensor basis is constructed on the following building blocks of Lorentz- and color-tensors. 
Color-tensors: 

^abcd ~ ^abScd, ^libcd ~ ^ac^bd, Cj^bcd ~ ^ad^bc, Cj^bcd ~ fabnfcdn, 
^ abed ~ facnfdbn 

Lorentz-tensors : 

From these, a preliminary generator system of the tensor space can be constructed 



Bi = 




B2 = 




^3 = 


= L<^'^C^3) 


Bi = 




B5 = 




Be- 


= i^^'qi) 


B7-- 


= L(^^C^2) 


B8 = 




B9 = 




Bio 




Bii 




B12 ~ 


- L(^^C(2) 


B13 - 




Bi4 




Bi5 
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with the Lorentz/color- indices left implicit. The tree- level tensor-structure of the four-gluon vertex 

is not a member of the preliminary system. To construct a system containing the trcc-lcvcl structure a Gram- Schmidt- 
algorithm is applied. It yields an orthogonal system of basis-tensors ^■'^T^hat- projection purposes it is also useful 
to define normalised quantities {j)U^^J^'^ such that 



%hcd 



y'' ^abbcdj ' -'-abed — abed ' ^ abi 

= Sjk. (Al) 



The basis system constructed this way is given by: 

^(1) ^ i08Arj(jVg-l) (~-^4 + 2B5 -I- 2Bg - BiQ - Bi4 - B15) 

U(2) = 468JV^(V^-1) (~-^4 + 2B5 -I- 2Bg — BiQ + 5Bi4 — B15) 

rr _ 1 

(3) 3510jVg(JVg-l)+ i''^™! (JVg-l)''-21060JVc(JVg-l) 

■{B4 — 2B5 — 2Bq + Bio + ^-Bis -|- ^Bi4 + B15) 

t^(4) = a7||M ( Ar2 _ 1)+ (iVf -1)2 - JVe( JV2 - 1) - M|p (;vj - 1) 

/ o or? O D I D I 306 o 693 r> 263 o 1 D ^ 
•^i>4 — /iiJs — ZlJg -h ±!io -|- —-D12 4o"-°13 lo"^14 + -°15j 

7-T _ i 

(5) 5616JVg (JVg - i"ia«2ua ( jyg - 1)2 _33696JVe (JVg - 1) - ^^^^^ ( JVg -1) 

•( — 8i?4 -|- I6B5 -|- I6B9 — 8i?io -|- — ^Bi2 — ^Bi3 -|- Bi4 — 8-B15) 

TT _ ^ 1 

'(<') 13816Ar^(iV^-l)-|-294(Ar^ -1)2 -1372 Ar„(iV^-l)-294(Ar^-l) 

■{B4 — 2Br, — 2B9 + ^Bio — ^Bii -|- jBi2 -|- I-B13 — ^-^14 + -B15) 

TT _ 1 

(7) 3455aii jv2(jV2_l)_|_9g(jV2_l)2_44gjV^(jV2_l)_96(Af2_l) 

/ 136 R 272 D I 2776 r> 104 d OR_lRj_R 1244 r j. 136 r \ 
■[ — -04 9--D5 + ~90~-°9 9~-°10 ~ ^-Dll + -D12 + -013 90~-°14 + —^15) 

U(s) = ( ivg ^"^^f ^'^ ( Afg - 1)^ - ^^^f^ iVc ( Afg - 1) - ^ (ivg - 1) 

/' R -L O R -L 846 R I 193 R I 178 R 18 R _L 9 r 81 r 68 r 1 r \ 
•(, — -04 + ^-05 + ~5~-D8 + ~5~-D9 + ~5~-DlO 5" -Oil + 5-°12 2"-°13 5"-°14 + -Disj 
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f/(9) - i.6U^i3'i jV2(JV2_l)+ iui'MMi (jV2_l)2_ jV^(Af2_l)_ 4«-i6^m4 (jv2_-Y) 

(' R _L R _i_ 11304 R 2844 n 5606 r i 178 r 

•(.-i>4 + ^-D5 + -5o--t>7 50" -^8 50~-^9 + ~^10 

18 r> 2736 r i 801 r , 1204 r n \ 

~T^11 ~ —^12 + -50"^13 + — -fc^W - -«15j 



'{10) - 4te4041ei jv2(jV2_i)+ I'^ag^^Oray (jv2_i)2_ 4elB4|Ma jv^(jv2_i)_ 4ij^^a^m4 (jv2_— 



('OR /I R _L '0<4 R laab R i38b R i iY8 R J5b R 



50 50 50 ' 50 

'^^^11 + 1^12 + 1^13 ~ -^14 + 2i?i5) 



no 



t^(ll) — 4aoati 



f^Af2(Af2-l) + 1440(Af2-l)2-6720Afc(A'^-l)-1440(7V2-l) 



/'107r i25rior R R 37 r 11 r 
■V-6~^4 + —-D5 + ^^6 — i37 — £>8 — 3Q-D9 

—8B11 + AB12 + 4i3i3 — -nr^i4 "T^^is) 



1 

56707V^(Af^-l) + 1440(Af^-l)^-6720Arc(iV^-l)-1440(Af^-l) 



•(^-64 ^ ^^5 ^ 2_B6 + + -Bs ^ + T^io 



+ 8-B11 — 4Si2 — 4i3i3 g--Si4 6~-^15) 

^(13) ~ 1080JV^(JV^-1)+7290(A'^-1)^ -6480^0(^^-1) 

•{^Bs + 5B4 + 5B5 — — Bg — Bio — 1^13 — Bi4 — B15) 

TT _ ^ 1 ^_ 

^(14) 2520Af^^(Ar^^-l) + 13666(7V^-l)^-14976Arc(A'c^-l)-6048(Ar^^-l) 

•(3OB2 - f S3 - I5B4 + 5S5 - 6B7 + 1^8 + 3^9 - Bio 
— QB12 + |i?13 + 3i?i4 — Bis) 

7-r _ 1 

^(15) ^ 1920Ar2(Ar2_i)+9720(Ar2-l)2-ii520Afc(A'?-l)-6480(Af2-l) 

•(25Bi - 5^2 - 5B3 + f B4 - f B5 - 5^6 + B7 + Bg - 1^9 + |B 



-bBii + i3i2 + -B13 — I-B14 + 1^15 
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APPENDIX B: THE CONSTRUCTION OF THE BOSE SYMMETRIC TENSOR BASIS 

To construct a basis of Bose-symmctric tcnsor-structurcs out of the structures given in appendix lA} one first has to 
construct a matrix representation of the permutation group with respect to the tensor structures. A general tensor 
in the Unear space V given by the tensor-basis from appendix [X] is represented by a vector 

15 

T = ^a,. «[/:^^/, 
1=1 

EE (ai,a2, . . . ,ai5) (Bl) 

An important feature of this tensor-basis, is that the tensor-structures constructed from it are closed under permu- 
tations of the external momenta. I.e. no tensor-structures not included from the very beginning are created by such 
permutations. This also means that permutations of the external momenta map the tensor space onto itself. Let V 
be a permutation of the external momenta, then 



■y e V : Vv eV. 
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This means that there is a matrix-representation of the permutation. For the four-gluon vertex, there are 4! permu- 
tations. Let M(^j^ be the matrix-representation for the j-th permutation. What one is interested in are the vectors 
V G V that are invariant under all permutations. To find them one first has to find the eigenvectors with respect to the 
eigenvalue one for each member of the matrix-representation of the permutations. The eigenvalue one in general can 
have different geometrical multiplicities for each of these matrices. This restricts the dimension of the Bose symmetric 
tensor space. The maximum number of linear independent vectors in the Bose symmetric space is the difference of the 
lowest and the second lowest geometrical multiplicity of the eigenvalue one in the set of matrices of the representation 
of the permutations. 

For the four-gluon vertex tensor-structures the lowest geometrical multiplicity of the eigenvalue one is five, while 
the second lowest is eight. Thus one ends up with an upper limit of three dimensions for the Bose symmetric tensor 
space. Having calculated the eigenvectors of the representation matrices, one can construct the full Bose symmetric 
linear tensor space. Let ek,i be the l-th eigenvector with respect to eigenvalue one of the fc-permutation matrix. 

Consider two permutations Vi and V2 , with matrix representations Mi and M2 ■ Denote the eigenvectors with respect 
to the eigenvalue one of these permutations as Vi and Wi respectively. Let the geometrical multiplicities of the 
eigenvalue one be /ii and /X2. The vectors, that are simultaneously included in the eigenspaces of two permutations, 
are found as the solution of the equation 

Y.^^■V^^Y.~^J■WJ. (B2) 

i=l i=l 

If this coupled system of algebraic equations is determined it can easily be solved by standard methods. This is the 
case, when there is only one eigenvector that is invariant under all permutations. The case of an under-determined 
system also is not problematic. The overdetermined system (which occurs in the case of the four-gluon tensors) is 
more complicated. It is then useful to reformulate the problem. Let A be the matrix, that is constructed from the 
column eigenvectors of both permutations in the following way 




With this definition 

A- (fi,..., e^,, C"^J^ = (B4) 

is equivalent to eq. (jB2p . The non-trivial kernel of the matrix A consists of -\- /Z2-dimensional vectors, whose first 
Hi components are the solution for the ^1 and the others are the solutions for the . The kernel of the matrix can 
then be evaluated using standard methods. For the four-gluon vertex the result is given in the main body of this 
work, eq. PH)) . 

APPENDIX C: DECOMPOSITIONS OF HIGHER GREEN'S FUNCTIONS 

In the following the decompositions of the reducible vertices in eq. into one-particle irreducible vertices will be 
given [2^, [3^. The reducible functions are denoted as T and the irreducible as F. 

1. The decomposition of the four gluon reducible Green's function: 
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APPENDIX D: ANALYTICAL INTEGRATION METHODS 
1. Collected integrals 

1. From ref. kal we use: 



(a) 



(b) 



sin^^'-^xdz f l\ f 1 1 



(1 + 2acosx + a^j \ 2/ \ 2 2 



[Re^ > 0, \a\ < 1] 



1 

(1 - x)'^^^x''^^ pFq{ai, . . . , ftp; 5i, . . . , 6,; ax)dx 







p+\Fqj^\(y, ai, . . . ,ap; + v, bi, . . . ,bq;a) 



r(/i + J/) 

[Re^ > 0,Rei/ > 0,p < g + l,ifp = g + 1, then \a\ < 1 



(Dl) 



(D2) 
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2. By partial integration, it can be seen: 



M a sin^'--^ 9 av f sin^^+^ 

dOCOaO X --77 = / dO rrr 

{l + 2acos9 + a^Y nJ (l + 2acoae + a'^Y+^ 





[Re/z > -1, \a\ < 1] 



(D3) 



2. Integrating with the Chebyshev-expansion 



A continuous function can be expanded in a series of polynomials. Most commonly one expands the function in 
a Taylor-series. But the convergence properties of a Taylor-series can fail to be sufficiently good, since the error 
of the approximation can be concentrated in a special region of the considered integration interval. When dealing 
with integrals on finite intervals, the Chebyshev-expansion can be an alternative, ft has the advantage that the 
approximation error is smeared out over the interval. The integral of the original function is reduced to integrals over 
Chebyshev-polynomials. 



The Chebyshev-polynomials are 

Tn{x) = cos(n arccos(a;)). 
A function f{x) can be expanded over these polynomials 

N-l 



Co 

2 ' 



(D4) 



(D5) 



with the coefficients 



N 



cos 



k=l 



j{k-l/2)7T 
N 



cos 



(fc- l/2)7r 
N 



(D6) 



are the zeros of the n-th Chebyshev polynomial. 



where N is the order of the expansion. The abscissas cos ^^^^ 
For a more detailed discussion see [46l |. 

To integrate a function using its Chebyshev-expansion, one transforms the variable, so that the integral is on the 
interval [-1,1]. 



dxf{x) 



x=r(y) dr 



dr r 

-1 dyY^c,TM-Y 

N-1 1 



dr 
dy 



-C0 + 



(D7) 



Thus the integration has been reduced to an integration over Chebyshev-polynomials. The integral over the 0-th 
Chebyshev-polynomial yields 2, while the integral over the 1st vanishes. 
The integral over the j-th Chebyshev polynomial yields 



cos(j7r) + 1 



dx Tj (x) = 

1-7 



(D8) 



Thus one obtains 



N-l 



cos(j7r) + 1 



Plugging in eq. (D.6) and denoting yi- = cos (^ ^''~]^^^^ ^ one finally gets 
5 

l^^fi^y- Ndy 



k=l 3=2 fe=l ^ ^ 



1 



3. The UV Integral 
In the UV analysis of the Ghost- and the Gluon-Box integral of the form 



oc 



(^iog(!) + i) 

yn 



occur. To evaluate them, one substitutes z = to log (f ) + 1 yielding 



I = 



us 



n—1 



(ti-l)z 

dze "2° 



z{_x) 



and then ^ 



(n-l)z 



. One gets 



e " 



az(x)) 



■. + l,(n-l) ■og(£) + 



s"-i(n- 1) 

Since a; is large in the region of interest, one can employ the asymptotic expansion 

TM-l 



r(a,a;) =a;«-ie-^ 



E^^-:# + 0(N"") 



,m=0 



with the Pochhammer-symbol (a)r 



(a)o = l, = ^^Y{a)^ = a{a + l) . . .{a + n-l). 

Keeping only the first order of the series, one results in 



/ 



dy- 



(n - l)a:"-i 
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